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and real part of all deeply virtual Compton scattering amplitudes. We note that a consistent 



(N 
> 

in 
o , 

[ We calculate spin, charge, and azimuthal asymmetries in deeply virtual Compton scattering 
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O at leading twist-two level. The measurement of these asymmetries gives access to the imaginary 
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Q^i description of this process requires taking into account twist-three contributions and we give then 
a model dependent estimate of these asymmetries. 
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1 Introduction. 



The scattering of electroweak probes off hadrons serves as a clean tool, free of complications 
from hadron-liadron reactions on both the theoretical and experimental sides, for the extraction 
of reliable information on the substructure of strongly interacting particles. Using the photon 
as a probe, (the absorptive part of) the forward virtual Compton (VC) process 'y*{q)N{P) 
'-f*{q)N{P) allows to study the strong interaction dynamics and at the same time it has a simple 
QCD description in the hard regime, — when — ^ '^Ldr- more general settings it is 
instructive to address the non-forward scattering 7*^*^(gi)A^(Pi) 'J^*\(l2) N {P2) . The systematic 
approach to its calculation is established only for the deep inelastic domain, where the QCD 
factorization theorems separate short and long distance phenomena into a perturbative parton 
subprocess and soft functions which encode information about the strongly coupled regime. The 
latter, known as the skewed parton distributions (SPDs), being studied for some time |^, have 
attracted increased attention in light of the conceivable opportunity to learn more about the 
spin structure of the nucleon 0. They are also of interest in their own right being hybrids of 
parton densities/distribution amplitudes and form factors. They share properties of the former 
in different regions of phase space which have been studied perturbatively to a great extent at 
one and two- loop orders 0]. Deeply virtual Compton scattering (DVCS), with = 0, proves to 
be an experimentally accessible reaction |^. In electroproduction processes of a real photon there 
is a strong contamination of DVCS by the Bethe-Heitler (BH) process. In view of the extreme 
interest to extract, or at least to constrain, SPDs it is timely to address the question of the best 
observables that allow to get rid of unwanted background. Fortunately, the interference of the 
two processes provides a rich source of information. It was suggested that diverse asymmetries 
P, H, 1^, H, P, |l^ can be used to disentangle the real and imaginary parts of DVCS and thus give 
access to SPDs. In the present contribution we consider a number of spin, azimuthal and charge 
asymmetries which share these properties and give predictions for the kinematics of the HERA 
and HERMES experiment. Spin asymmetries make it possible to extract the imaginary part of the 
DVCS amplitude and thus, due to the reality of SPDs, which holds owing to the spatial and time 
reversal invariance of strong interactions, give directly a measurement of the shape (at leading 
order in in complete analogy to DIS) of SPDs on the diagonal t = ^. 

The paper is organized as follows. In section ^we calculate the squared amplitude for DVCS, 
BH, and the interference term to leading twist-two accuracy. In section |^ we present simplified 
formulae for different cross sections and give a qualitative discussion of the feasibility to measure 
the DVCS leading twist-two amplitudes in different kinematical regions. After introducing models 
for SPDs, we give an estimate for different asymmetries for HERA and HERMES kinematics in 
^. Finally, in section ^ we summarize. 
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(a) (b) 
Figure 1: The virtual Compton scattering amplitude and the Bethe-Heitler process. 



2 Cross section. 

To start let us discuss different contribution to the differential cross section of the electropro duct ion 
process e{k, \)N{Pi, Si) e{k', \)N{P2, S'2)7(g2, A), given by the standard formula 

The scattering amplitude squared |Tp in the cross section, contains beside the VCS [Fig. 0(a)] 
and BH [Fig. ^(b) and crossed contribution] parts also the interference term: 

A',52,A 

The BH-amplitude is purely real and is given as a contraction of the leptonic tensor, at leading 
order in the fine structure constant a, 

V = Hk', A') [7m(^ + lu{^' A), (3) 

with the hadronic current 

J, = t/(P2,^2){Fi(A2)7. + ^F2(A2)fx..^|t/(Pi,^i), where A = P2 - Pi = gi - ^2, (4) 

parametrized in terms of Dirac, Pi, and Pauli, P2, form factors normalized according to Pf (0) = 1, 
Pf (0) = Kp = 1.79, and Pi"(0) = 0, P2"(0) = k„ = —1.91, for proton and neutron, respectively. 
Thus, the BH amplitude is of the form 

Tbh = -^.elL^^Ju. (5) 

The form factors are known fairly well from experimental measurements and can be parametrized 
by dipole formulae in the small region 

GU^') = (1 + Kp)-'Gl,{^') = n-'Gl,iA') = f 1 - ^) , ^^(A^) = 0, (6) 



where we have introduced the electric, G'^(A^) = ^{(A^) + ^^F2(A^), and magnetic, G\^(A^) = 
FKA"^) + F2(A^), form factor characterized by cutoff mass my = 0.84 GeV, see e.g. |jl3|. The 
hadronic tensor Tycs is: 

p'^ I — for p~ 

Tvcs = T-e^.T^'uik'Mk), where , , , • (7) 
qf ^ [+ for e+ 

It is defined by the time ordered product of two electromagnetic currents 

T,,{q,P,,P,)=z J dxe'^-''{P2, S2\Tj,{x/2)U-x/2)\P,, Si), (8) 

where q = {qi + Q'2)/2 (and the index /i refers to the outgoing real photon). It contains for a 
spin-1/2 target twelve^ independent kinematical structures 0. In this paper we restrict ourselves 
to the twist-2 part of T^i, that does not contain transversal photon spin flip contributions. Such 
contributions arise in the next-to-leading order of perturbation theory due to the gluon transversity 
^ [Tl| and are seperately considered in . From the structure of the OPE we immediately learn 
that these contributions are contained in the following form factor decomposition^: 

T^uiq, P, A) = -^^^^^ - ie^^g^^ + ■ ■ ■ , (9) 

P.q P.q 

where P = Pi + P2 and the gauge invariant tensors t^,^ = V^ptp^Vau are constructed by means of 
the projection tensor V^^u = g^i, — qip,q2u/ <li-<i2- The ellipsis indicate twist-three and higher contri- 
butions. The vectors Vi^ and axial-vectors Aiy can be expressed by a form factor decomposition 

Vi, = U{P2, S2) [Hi^, + £1^^) U{Pi, 5i) + ■ ■ ■ , (10) 

Ai, = U{P2, S2) (^i7m75 + £1^) UiPi, 5i) + ■ ■ ■ , (11) 

where higher twist contributions are neglected. These form factors depend on the following vari- 
ables 

P.q 4 

Note that in general (for off-shell final photons) a second scaling variable 77 = ^ appears, which 
is related however to ^, i.e. r] = (l — ^ — ^. The amplitudes are given as convolution in 
t, (g) = / dt, of perturbatively calculable hard scattering parts with SPDs: 



} (e, g^ A^) = Ti(e, g^ /^', ® { f } (t, e, A^ /x^), (12) 

I (e, Q^ A^) = fi(e, Q^ /i^ t) ® 1 1 1 {t, a^, f,^), (13) 



^12 = i X 3 (virtual photon) x 2 (final photon) x 2 (initial nuclcon) x 2 (final nuclcon). The reduction factor 
f/2 is a result of parity invariance. 

^We adopt throughout the conventions of ||lj], e.g. e°"^^^ — 1 
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Figure 2: The kinematics of the reaction e{k)N{M) — ^ e{k')N [P 2)'^{q2) in the rest frame of the 
target. 

with summation over the different parton species imphed and yU^ being the factorization scale. 
The hard scattering amphtudes are available in next-to-leading order (NLO) approximation and 
they read in LO for a quark of charge Qi 



^j.i(0) ^^^^^ 



(14) 



1-t/^- ie 
with — (+) for parity even (odd) cases. 

In the consequent presentation we give our results in the laboratory frame (see Fig. ^) in which 
we use the kinematical variables 

k = {E,0,0,E), A;' = (E',E'cos0eSin^e,E'sin0eSin^e,E'cos^e), (15) 
Pi = (M, 0,0,0), P2 = {E2, IP2I cos0Ar sin^iv, IP2I sin 0Ar sin IP2I cos 6'Ar). 

Furthermore, we introduce the azimuthal angle (pr = 0iv — 0e between the lepton and hadron 
scattering planes as well as (j)s = (pN + <Pe- The spin vector of the spin-1/2 target for longitudinal 
and transverse polarization is given by 



5= (0,0,0, A) with A = ±l, 5 = (O,cos^,sin^,0), 



(16) 



respectively. 

From the experimental point of view one works with the variables = —qf and x = 
—ql/{2Pi.qi), which are related to our variables by 



I J A2\ 1 „ 



X 



(1 



A2 
2^ 



X 



2-x 1 + 



A2 



2 - X 



(17) 



After performing the phase space integration we obtain for the differential cross section (|1|) 



da 



y 



da 



dxdQ^d\A^\d(l)r dxdyd\A^\d(l)r 8nQ* 



a^xy^ ( 4MV\"'/' 



^3 



(18) 



where we introduced as well the conventional variable y = Pi.qi/Pi.k [y = 1 — E' /E in the frame 

Here we present simple analytical expressions for the amplitudes entering the cross section for 
positron beam. Changing to electrons will generate a minus sign for the interference terms. To 
deduce them we present at first the amplitudes squared in terms of the form factors. The result 
for the DVCS amplitude reads in leading twist-2 approximation (we set |e| = 1) 

A(2_y) ^ (g.A q.V^ + q.V, q.A\) , (19) 

where A = 1 means that the spin of the lepton is parallel to the beam direction. It is obvious that 
the leading contribution scales as 1/ Q^. 

The exact squared amplitude for BH reads in terms of the electromagnetic currents: 

8 q.J k.J^ + k.J q.J^ - q.J q.J^ - 2k. J k.J^ 



BH 



A2 (2A;.A- A2)(Q2 + 2A;.A) 

~ A^ y + 2(2A;.A - A2)(Q2 + 2A;.A) ) 



4iA 
+ 



Q2 + 4A;.A-A2)e,Ajjt+A2 



A4(2A;.A - A2)(Q2 + 2A;.A) 
with the latter being defined in Eq. @). Instead of the exact expression we may use an approx- 
imated one. To be consistent we have to expand the squared BH-amplitude, which starts with 
1/A2, and the interference term up to the same order as the squared DVCS-amplitude, namely to 
1/ accuracy. However, one has to take into account that the lepton propagator in the w-channel 
of the BH amplitude gives a contribution which behaves as {k — ^2)^ = — (1 + 0{1/Q)). 

Therefore, the Taylor expansion is not legitimate for large y and sets, otherwise, an upper limit 
for y, namely 1 — y ^ M"^ j Q2. To avoid this problem, we have to expand the propagator. 



02 r _A2 A2 . a2 

{k-q^f = -^\l-y + 2J-—Jl-^JT^VT^cos'- 



y [ ^ V QM A2 V «^ ^^^.vY-r; 

+ {1-X)[l- 2^] -(l-y)L + 2(1 - x)^] cos(20,, 
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(21) 



and the remaining parts separately in Taylor series which results in a Pade-like approximation for 
the squared BH term. Here A^j^ is the minimal value of A2 which is defined by the kinematical 
restriction A^j^^ = —M'^x'^/ (1 — x + xM"^/ Q'^). Taking only the first terms in the expansion, thus, 
neglecting all contributions formally suppressed by l/Q, we get 

2 2(2-2y + y2) fjji ^ q.Jq.J^ \ ^ \{2 - y)y' te.^jj, 

- (l-y) I A^ + ^^^J~^ (l-y) A4Q2' (22) 
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which, however, is vahd for y <^ 1. 

For interference terms we adhere to the same approximation and give here only the leading 
contributions in 1/Q, i.e. 



^H^vcs 



1-y A2Q4 

A (2-y)i/ e 
1-y 



-^"r^A^ii^^-^^^j (^. + 2A.|^j,4-2..,..^J. (23) 

Note that for a consistent expansion up to order l/Q^ also twist-three contributions having both 
kinematical and dynamical origins must be considered (see Ref. |T5| for the case of a scalar target). 
In general these contributions are poorly understood and have to be studied in more detail. 

As a next step for the calculation of the cross section we have to sum over the spin of the 
outgoing proton and write the final answer in terms of SPDs and electromagnetic form factors. 
Here we give the results for polarized spin-1/2 target, where A = 1 means polarization along the 
lepton beam, averaged over the azimuthal angle 0^. The DVCS amplitude squared I'TdvcsP = 
+ I'^DVCslpob with pol = {LP, TP}, consists of the elements 
2{2-2y + y') 



2iek, 



kqAJ' 



I'^VCslunp 



IT |2 
I^DVCSlunp 



?/2(2-x)2Q2 



V + (2-x)^^)£,£r 



4(1 - x) (HiHI + HiHl) - x^ [HiSl + SiHl + HiSl + S^Hl 

o A2 



X 



4M2 



(24) 



ITr 



DVCS I LP 



DVCS I TP 



2AA(2 - y) 
\/(2-x)2Q2 



4(1 - x) [HiHl + HiHl) - x' [Hi£l + EiHl + HiEl + EiHl 
A2 



(25) 



8cos(<P-0^/2)(2-2y + ?/2)^/r^ /_A2 / 
7ry2(2-a;)2Q2 ^^PV 

< {(2 - x) [(Re7^l)(Im^l) - (Im7^i)(Re^i)] - x 

2Asin(^-0^/2)(2-y)v/r^ / A2 
7r?/(2-x)2Q2 V~M2 

-2(2 - x) (n,£*, + ni£i) + x^ (£i£l + £l£i 



A2. 

mm 



1 - 



A 



A 



(Re7^l)(Im£l) - (Im7^i)(Re^i)]} 
2x (ni£l + 'Hl£i) (26) 



For the BH term we use an analogous decomposition with 



I^H 



2 

unp 



IT |2 

KBhIlp 



IT |2 

I -'bhItp 



2(2-2y + y2 
(l-y)A2 
A\k{2-y)y 



1 



X 



x^ 



A2. ^ 

mm 



Fl + 2(Fi + F^f + 



A2 



A2. 

"mm 



1 F| 



[l-yW 



{Fi + F2) 



X 



X 



A2. 

mir 

A2 



Fi + Fi + F2 



(27) 
(28) 



8A sin(<f - (t)r/2) (2 - t/)y v^T^ ^ 
7r(l -?/)xM(-A2)3/2 ' 



A2 



A2 



(Fi + F2) (2xM2Fi + A^Fa) . (29) 
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For a polarized lepton beam the interference term is decomposed into a contribution for an 
unpolarized and an additional one for a polarized target, i.e. 2^ = ^h^vcs + ^vcs^h — 
^np(A)+2^oi(A), where 



^-^^^^ ^ ~ Vl-yyxV-A'^Q^ - ^cos(0.)Re|F,7i, + ^(F, + F,)7^, 

- m^^^\^M^r)^F.n. (30) 



4M2 'I yr^xV-A2Q2 



xRe{^(F, + F.) + ^.,) + - ^ (^F, + A (31) 



8(2-2|/ + |/2) 



cos (^ - 3(/)^/2) / A 



2 



(l-^)H 2F2(7ii + ?ii) 



27r V A2 / ' 'I 



-[(2 - x)Fi - xF^\E^ - xFa} ^ ^o^(^ + '/'^/2) + F,)Im{4(l - x) 

X (Hi - TYi) - - £i)} + - Im{2(l - a;)F2(7ii - H^) - [(2 - x)F, 

sin (^ - 30^/2) 



+a;F2]^i + a;(Fi+a;F2)£i} 
Ai. 



^ 8A(2-|/) 



X 1 



v/r^(2 - x)xV'Q2]^^ 
^ (1 - a;)Re{2F2(7ii + H^) - [(2 - x)F^ - xF^\E^ - xF^E^ 

'^^^ (^(F, + F2)Re{4(l - .) - H.) - x\E. - ~E.)\ 

j Rc{2(l - x)F2 {Ui - Hi) - [(2 - x)Fi + xF2]Ei 



(32) 



/ A^. 

I I mm 



+x(Fi + xF2)£i} 

3 Extraction of leading twist-two amplitudes. 

A strong motivation for the measurement of DVCS arises from the fact that the second moment 
of SPDs in the parity even sector is related to form factors appearing in the decomposition of the 
symmetric QCD energy-momentum tensor 0^^^. A gauge- invariant decomposition of the matrix 
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element of the angular-momentum tensor, M.^,y^a = x^Oa,y — x,yOa^, provides, therefore, a separate 
estimation of the total angular-momentum fraction carried by quarks and gluons 0. To achieve 
this goal it is necessary to interpolate the corresponding moments of the SPDs to forward kine- 
matics A = 0. The new information that is required (and not available from DIS) is contained in 
the SPDs appearing in the spin-flip amplitude Si. Having the spin puzzle in mind, we give special 
attention to the problem of extracting Si from measurements in different kinematical domains. 

The simplified explicit expressions ( ^^ - |3^ ) for the amplitudes squared allow us to discuss the 
extraction of the leading twist-two amplitudes from future experimental data. Moreover, they 
allow us to give a qualitative discussion of the ratio for the different cross sections in more detail. 
Note, however, that the formulae, presented below, are only valid in a small kinematical window 
and in general not sufficient for numerical estimates. In the following we divide the kinematical 
region into ^ A^j^, |A^jjJ < |A^| < M^, and |A^i^| < < |A^|. Furthermore, we separately 
consider the small x region. 



3.1 Small X region. 

Let us start with the small x region, which is for instance relevant for the HERA experiments 
We assume that the small x behavior of Si and Si is not one power less than that of Tii and Tii, 
and, furthermore, that the SPDs in the small x region are essentially determined by the usual 
parton densities at least in the DGLAP region. It is easy to establish that for small values of x 
the ratio of the real to imaginary part of the unpolarized amplitude Tii is about 0.3 (by means of 
dispersion relation the same magnitude has been obtained in Ref. H]). In the polarized case we 



find for instance for the Gehrman-Stirling parametrization that this ratio for TYi is of order 
1 — 1.7. However, the contribution of the latter is quite small as compared to Tii. 

Taking into account these numbers, we find that for longitudinally polarized target Eqs. 
can be approximated by: 

2{2-2y + y^) 

y2Q2 V X, ^^^2 



ITr 



DVCSI 



(33) 



_ 4AA(2 y) /j^-^^j^-^^ ^ ReHiReHi 
yQ^ ^ 



Contributions containing Si are down by two powers of x or they are proportional to x^2- The 
transversally polarized part offers an opportunity to measure Si. 



IT |2 ^ 4cos(^-(/.,/2)(2-2y + y^) -A2 Al^ vt P ^ a vr p M 

I^dvcsItp ~ V^TT^yl - -r^ [(ReHi)(Im£i) - (ImHi)(Re£i)] 



7r?/2g2 V M2 V A 



4Asin(^-0,/2)(2-?/) / A^ / A 



2 



1 - 



mm 



7r?/Q2 V M2 V A2 



(Im7ii)(Im£i) + {Reni){ReS 



(34) 
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The interference term starts with x~^, however, since the sea quark and gluonic contributions to 
the DVCS amphtudes are expected to grow with x~^, the DVCS cross section and the interference 
term have the same x dependence. For the longitudinally polarized case it reads for small x: 



sin(0r)Wl ^ f / \ 



X 



8cos(0,.)\/l - r 
X ^^(F, + + F,H,) } - _=X__Re|(2 - 2y + (35) 

The interference term might give in future an opportunity to access the imaginary and real part 
of the hnear combination FiHi — ^^-F2£^i due to the charge and single lepton spin asymme- 
tries for — A^jjj < — A^. For a polarized proton beam one may extract in an analogous way 
Im [x{Fi + ^2)7^1/2 + FiHij and in combination with a polarized lepton beam one also gets the 
real part of this expression. Note that one has now the opportunity, at least in principle, to extract 
the DVCS cross section (^) for unpolarized or double spin flip experiments and thus separate Tii 
from the £i contributions for the imaginary and real part. 

The squared of the BH term will generally start with (however, the spin dependent part 
goes only like x~^) and has therefore a similar x dependence as the other ones, i.e. xlmTii/Fi is 
of order one or so, and grows only slightly with increasing x. However, in comparison with the 
squared DVCS term it is multiplied by ?/^Q^/A^. Thus, one expects that the BH background is not 
large in the small y and — A^/ region. We conclude from our discussion that the measurement 
of the unpolarized and longitudinal polarized cross sections as well as single spin asymmetries it 
is feasible to disentangle the imaginary and real part of T^i, Tii, and £i at small x and — A^ of the 
order of one GeV^ or larger. For smaller values £i is in general kinematically reduced. Since £i is 
suppressed at least by one power in x, we conclude that this amplitude will not be accessible in 
the small x region. 

As mentioned, if — A^ starts to be smaller, the contributions proportional to £i die out. In the 
case of its value at the kinematical boundary, i.e. A^ ^ ^min^ see that Eqs. ( |27| - p9D are quite 
simple, and, for instance, for a longitudinally polarized target we obtain 

i{2--iy + v') . 4AA(2-y)y , 

Moreover, in this case the interference term for unpolarized or longitudinally polarized target 
drops out, too. It is worth to have a closer look at the ratio of the DVCS and BH cross section. 
The unpolarized cross section is essentially governed by the imaginary part of unpolarized parton 
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distributions: 

c^(^DVcs(A = 0, A = 0) (1 - y)x'^M'^ ImHl 



(37) 



It is remarkable that the ratio of double spin flip (DF) cross sections is proportional to the product 
of unpolarized and polarized parton distributions: 

^^DVcs(A = 1, A = 1) - d(7Dvcs(A = -1, A = 1) _ (1 - y)x^M^ ImT^iImT^i + ReHiReHi 
daBH(A = l,A=l)-d(7BH(A = -l,A = l) ~ y2g2 (Fi + Fs)^ ' ^ > 

For the sake of completeness we want also to mention that the interference term for a transver- 
sally polarized target does not vanish for ^min' while both the DVCS and BH cross section 

are reduced to the unpolarized ones. This interference term is proportional to 7ii — TCi. 

Note that at the kinematical boundary the spin flip contributions are generally suppressed by 
at least a factor x^. Thus, in this kinematical range it is hopeless to extract £i or Ei from any 
experiment. 



3.2 Asymmetries for |A^| > \^min\ 



Away from the kinematical boundary |A | 3> |A^i^| we can consider different cross sections to 
get separately information on the real and imaginary parts of the leading twist-two amplitudes, 
Til, . . . ,Si. Since the BH cross section depends only on the product of the Icpton and target 
polarization, i.e. on AA or A/((?) with f[$+7r) = —f{'P), we have a direct access to the interference 
term in polarized experiments to leading order in Moreover, the relative sign of the 

interference term is determined by the charge of the lepton beam. Thus, the following cross sections 
allow one to get access to the leading twist-two structure functions in the DVCS amphtude. In 
the approximation used above and neglecting terms O (A^j^/A^), they read: 

1. Polarized positron beam and unpolarized target: 

AsL^fJ = da^ - da^ (40) 
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2. Unpolarized positron beam and longitudinally polarized target: 



AsLNflfJ = dcr^ - da:^ = , . sin(0r 



(41) 



3. Unpolarized positron beam and transversally polarized target: ($ = {0,7r}): 



X 



cos(30,/2) ^^ _ ^^j^ + 7?,) - [(2 - x)F^ - xF,]S^ - xF,S^} 

-^^^^^lm{2il - x)F, {h, - H,) - [(2 - x)F, + xF,]8, 

+x(Fi + xF2)^i} 



(42) 



dM. 



4. Charge asymmetry in unpolarized experiment: 

unp j-_unp_ lQ{2-2y + y'^)^Jl- X 



a (T = 



./^^yx^/-^^Q^ 



COs(0r) 



(43) 



xRe <! FiTii + ;^(Fi + F^)^,^ - j^F^Si } dM. 



X 



5. Charge asymmetry in double spin- flip experiments with longitudinally polarized target: 



Ag^^da = d+a]-d-a\-Acda''-P^^^^J^^^cos{<f>r) (44) 



xRe{^(F, + F.) (n. + |.,) + F.n. - ^ + ^F.) .~ j.A.. 

6. Charge asymmetry in double spin- flip experiments with transversally polarized target: 

8(2 - y) 



A^^^da 



X 



d+ai - d'at - Acci(j""P = , , 

^/^^(2 - x)a;x/Q2M2 

sm(3(/.,/2) ^^^ _ ^)j^g|2^^(^^ + _ - x)Fi - xFaj^i - xF^Ei} 



27r 
sin (0r/2) 



TT 



Re{2(l - a;)F2 (h^ - Hi) - [(2 - x)F^ + xFs]^! 



(45) 



+x(Fi + xF2)£:i} 
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where dA4 = (l + ^^^^f-) ^ dxdyd\A'^\d(l)r. Note that for a consequent numerical treatment 
the A^in/A^ dependence cannot be droped and can be easily restored from Eqs. (|30H32D. 

As we see from this list, the two single spin cross sections of longitudinally polarized beam or 
target give us information on the imaginary part of two linear combinations of the four leading 
twist-two amplitudes. The real part of these two quantities is accessible via charge asymmetry in 
unpolarized or longitudinally double spin flip experiments. Obviously, at low |A^|, i.e. compared 
to 4M2 ^ 4 GeV^ there is a suppression of contributions proportional to E and E SPDs, which 
are theoretically not well constraint: the contribution proportional to E SPDs completely drops, 
while E is suppressed by x^. For the kinematics of present experiments — A^ ~ IGeV^, and thus 
it is not a good approximation to drop afore mentioned terms, A^/4M^ ~ However, in 

certain asymmetries these contributions are accompanied by a factor x and thus can be safely 
discarded at smaller values of x. Further constraints separately for the imaginary and real part 
arise from spin and charge asymmetries for a transversally polarized target. Fortunately, due to 
different angular dependences, we have indeed two further independent linear combinations of 
the DVCS amplitudes. But the kinematical pre-factors and an additional angular dependence 
of the denominator for large ?/, which we have dropped for simplicity, makes their practical use 
questionable. 

Finally, the whole cross section for unpolarized and polarized 7 production tests the leading 
twist-two approximation used above. However, to be consistent one has then to expand the BH 
and interference terms up to order 1/ Q^. As outlined in section § this is straightforward for the 
BH cross section, while in the interference term new twist-three contributions will enter. These 
additional terms will be worked out in a forthcoming paper. Here we should only mention that 
the different azimuthal angle dependence can be used to pick up different combinations of helicity 
amplitudes, entering at different twist levels. For instance, for small y it is justified to use the 
derived formulae which imply that all single spin and charge asymmetry cross sections integrated 
over the azimuthal angle 0^ vanish at twist-two level, for instance, 

4 Numerical estimates. 

To give phenomenological predictions for the asymmetries discussed above we have to specify 
models for the SPDs. Definitely, this is the main source of uncertainty for the numerical estimates 
we present in this section. (Of course, the primary goal of experiments is rather to constrain the 
skewed functions via the theoretical formulae.) 

Let us discuss in turn spin non-flip and flip functions. In the former case SPDs have a definite 
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limit for forward kinematics when they reduce to the famihar parton densities. For spin non-flip 
SPDs we choose an oversimplified factorized form of the and (t,^) dependence of the SPDs: 



W{t, e, A^, Q2) = FliA')q\t, e, Q'), H\t, ^, A^ Q') = G\iA')Aq'{t, ^, Q'), (47) 

where i denotes the quark flavour. Here F{(A^) and G'*(A^) are elastic parton form factors 
normalized to unity at the origin, -F^*(0) = G\0) = 1, and qit, ^, Q^) as well as Aq{t, Q^) are the 
non-forward functions specified below. The support of these functions is [—1, 1], where for t > we 
have the quark distribution and for t < the antiquark distribution, i.e. q(t, ^, Q^) = —q{—t, ^, Q^) 
and Ag(t,e, Q^) = Aq{-t,^, Q^). The normalization of H%t,^,A^, Q^) and H\t,i,A^, Q^) at 
A2 = 0, ^ = is determined by parton densities H%t, 0, 0) = g'(t), H%t, 0, 0) = Aq\t). Since the 
^ dependence drops out in the first moment we can constrain the A^ dependence by sum rules, 
e.g. 

f^dtH\t,i,A') = Fl{A^), f^dtE\t,^,A')=F^{A'), (48) 

with the Dirac and Pauli form factor, respectively. For non-polarized SPDs the valence u and d 
quark form factors in the proton can be easily deduced from (^ via " = 2 ^ + (^^"^^ Ff 
which results in 

2F^{A^) = 2Ff{A^) + Ff{A^), Ff{A^) = Ff{A^) + 2Ff{A^), for J = 1, 2. (49) 

For s (or in general, sea) quark contribution in the parity even sector the first moment vanishes 
and the sum rule ( ^Sf ) does not give any constraint. Nevertheless, the counting rule for elastic 
form factors tells us that for large A^ we have a (A^)^^ behaviour. This suggests the following 
dipole fit with the mass cutoff my chosen as for valence quarks: 

(A^) = — ^G1T(A2) = f 1 - ^) ' , (50) 

where Kgea will be specified below. 

So far the model has been governed by the known forward densities. Unfortunately, a similar 
reduction is absent for the helicity flip amplitudes. For E^ we adopt nevertheless 

E\t,i,A^)=r\t,i)r^{A^) with r\t,i) = q\t,i). (51) 

The identification of r^{t,C,) with q'^{t,C,) ensures the sum rule for valence quark contributions. 
Note that the parameter Ksea, normalizes the sea quark contribution, for instance, Kgea. = provides 
^sea _ g j^-^ ^^iq axial vcctor channel the quark form factors can be read off from the iso-triplet 
axial form factor G'P''(A^) of the /?-decay and related by isotopic symmetry to the form factor in 
question. The decay constant g^2^ is expressed by Goldberger-Treiman relation ^f^'' ^ -^^^fngnNN 
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in terms of the pion decay constant /^r, nucleon mass M and vrA^A^-coupling, g^NN, and has 
the numerical value g^^^ = 1.267 DTS). If we assume the same A^-dependence for the iso-singlet 
^{^^(A^) with the same cutoff and a constant g^^^ we get for quarks 

Gi(A')= (^1-^^ , for z = K,4}, and GHA^) = (^1 - (52) 
with the scale = 0.9 GeV pi. 



Finally for the polarized spin-flip amplitude it was observed [T^, [18], |T^ that, similar to the 
/9-decay effective pseudoscalar form factors [|T^, one can approximate E at small A^ by the pion 
pole so that one ends up with the model 

1 4r7?'*A/f2 4 

= E'' = -e{^-\t\)<PAt/09A^'), with g^{A')= ^ and = -(1 - x^). 

24 — A^ 6 

(53) 

It is reliable to assume that the SPDs in the DGLAP region can be modelled by the forward 
parton distributions measured in inclusive reactions. In the simplest case we assume that this is 
the only contribution. We refer to this model as the forward parton distribution (FPD) model 
which has no skewedeness dependence at the input scale Qq: 

{A)q\t,tQl) = {A)q\t,Ql) for t > 0. (54) 

The contributions for t < are easily restored by means of symmetry. Although the input 



q{t, Q') = /' dx /' dy6{x + ^y- t)f{y, x, Q^). (55) 



distribution does not depend on ,^ a small evolution step does generate such a ^-dependence 
A further model is based on an proposal for the so-called double distribution (DD) |]l|, namely, 

The functional dependence in the x-subspace is given by the shape of the forward parton density, 
f{x), while the y/[l — |x|]-dependence of the integrand has to be similar to that of the distribution 
amplitude. Thus, f{y, x) is given |]TB[ by the product of a forward distribution f{z) (more precisely 
q{z) for quarks and zg{z) for gluons) with a profile function n 

f{y,x) =7r{y,x)f{x), (56) 

where vr for quarks and gluons is given by 

- 4 [1- NP ' ^ - 16 [1- W • ^''^ 

Now we are in a position to give numerical estimates for the cross sections defined in the preceding 
sections. 
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4.1 HERA kinematics. 

In the small x kinematics, 10~^ x < 10^^, first estimates for the unpolarized azimuthal and 
single spin asymmetry has been given in the factorization p[ and BFKL approach. In the 
following we evaluate numerically different asymmetries in order to demonstrate that Tii, Tii 
and Si are measurable in HERA experiments. In the following we deal with the FPD model, 
where we equate the SPDs with the parton densities taking the MRS A' and the GS A pi 



parametrization at the input scale = 4 GeV^ and using u = d = s/2 as well as Au = Ad = As. 
For simplicity we do not discuss the dependence of our predictions due to the perturbative 
evolution. 

The unpolarized azimuthal asymmetry is defined by 



and its first signature has been seen by the ZEUS collaboration [^. In the approximation ( PT]) the 
unpolarized squared BH term is ipr independent. However, it turned out that this approximation 
may provide misleading results due to the neglected azimuthal dependence of the BH process. 
In Fig. |^(a) it can be seen that the dropped terms cause a strong 0,.-dependence and that only 
for the Pade-type approximation introduced in section ^ [see Eq. ( |2T| ) for the expansion of the 
propagator] there is a good agreement with the exact expression. Since the BH cross section is in 
general suppressed in the upper hemisphere, one may expect quite different predictions depending 
on the approximations involved. Indeed, for small values of x, where y = Q'^/xs with the center 
of mass energy s = 4 ■ 27.5 ■ 820 GeV^, we find a strong deviation, even for small — A^. Note that 



the interference term is only taken into account up to order 1/a/— A^Q^ and at present it is not 
clear whether the 1/Q^ term is crucial or not. 

To get information on |-Fi7Yi — ^^-p2^^i} in a cleaner way one may use the azimuthal asym- 
metry of the charge asymmetry A'^'^da = d~a'^^^^ — (i+cr"°^ defined in Eq. ( ^31) (in the following 
we restore the A^j^/A^ dependence in the interference amplitudes): 

and the single (lepton) spin asymmetry with unpolarized target Asl = da^ — da^ defined in Eq. 
iOl): 



AsL = T ^^/i ■ (60) 

The former (later) one is proportional to the real (imaginary) part of the linear combination that 
was mentioned. Both of them give sizeable effects. Although is proportional to the imaginary 
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Figure 3: (a) Ratios of different approximations of tlie BH cross section to the exact one in depen- 
dence of the azimuthal angle for = 4 GeV^ and x = 10~^ for the Pade-type approximation 
for = —0.1 GeV^ (sohd) and = —0.5 GeV^ (dashed) and for the leading approximation 
( p7D (dash-dotted and dotted, respectively), (b) Unpolarized azimuthal angle asymmetry for 
A^ = — 0.05(— 0.25) GeV^ versus x. Solid (dash-dotted) and dashed (dotted) lines show the result 
for the Pade and leading approximation, respectively, where Kg^a is set to zero. 



part which is growing as x~^ this asymmetry is suppressed by a kinematical factor y as compared 
to the charge asymmetry. Therefore, we find in the considered kinematics a two times larger value 
of the charge asymmetry as for the single spin one. This should be reversed for larger values of y 
(see also the discussion in Ref. |]10|]). In Fig. ^ we demonstrate the influence of the amplitude Si 
which is quite sizeable for larger values of — A^, where — A^/ remains small. 

To measure Tii one have to consider the proton single spin asymmetry Asln = dc^i — d<^i given 
in Eq. (|l]): 



A 



■J da 



-I 



2n 



SLN 



d<f)r 



27r 



(61) 



and the charge asymmetry in double spin-flip experiments, i.e. Aq da = d <t^ — d~^a^ — Acda^'^^ 
defined in Eq. (|3) 



4DFL 



aDFL j„ 



Jn 



2tt 



(62) 



Unfortunately, for the considered value we find for our model that this single spin asymmetry 
is compatible with zero. Although, we have an enhancement for small y due do the factor l/y in 
comparison to the lepton single spin asymmetry, this enhancement cannot compensate the weaker 
rise of Tii for small x. For larger value of Q, i.e. also larger x, and — A^ and small y we found a 
few percent effect. Note also that the ratio of proton to lepton single spin asymmetry gives us the 
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Figure 4: Unpolarized azimuthal charge asymmetry Aq for = 4GeV^ and x = 5 ■ 10~^ 
with Ksea = "~2 at LO for the complete expression (sohd curve) and neglecting the E\ contribution 
(dashed curve) plotted in (a) versus — A^. The same asymmetry for = —0.05 GeV^ (solid curve) 
and A^ = -0.25 GeV^ (dashed curve) at LO is shown in (b) for the region 1-10 ^<x<2-10 ^ . 
In (c) and (d) the same is shown for the electron single spin asymmetry. 

ratio (Im7ii)/(Im7ii) times l/y, which is sizeable at not too small x. The azimuthal asymmetry 
in double spin flip experiments normalized to the unpolarized cross section is tiny. However, if 
we would compare it with the extracted double spin flip part, we would get quite sizeable effects 
which are of the order of 10% or even more. Furthermore, to explore ?ii separately, one can make 
use of the discussion given in section p.l| . Especially, for small y the BH cross section is suppressed 
by ip' in comparison to the DVCS one and, thus, the former one drops out. 

Let us note that the perturbative NLO corrections to the imaginary part for small x are of 
order of 20% for each quark species as well as for the gluon contribution. 



4.2 HERMES kinematics. 

Now let us turn to the HERMES experiment with a = 27.5 GeV positron beam scattered 
by a hydrogen target and give predictions for the asymmetries which can be accesed there. To 
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Figure 5: Perturbative leading order results for the charge asymmetry for an unpolarized beam 
(a), single spin asymmetries for a polarized positron beam (b) and an unpolarized target; as well 
as for an unpolarized lepton beam and a longitudinally (c) (transversally (d)) polarized proton 
target versus x, for = 4 GeV^. The predictions for the model specified in the text are shown 
as solid (dotted) curves for A2 = -0.1(0.5) GeV^ respectively. The same model however with 
neglected spin-flip contributions are presented as dashed (dash-dotted) line for the same values of 

give numerical predictions we take the model for SPDs deduced from the double distribution 
model and forward (un-) polarized parton densities from Ref. (p3[) with Kgea, = —2. As a 
starting point we choose = 4 GeV^, x-range 0.1 — 0.4, and the t-channel momentum transfer 
— = 0.1 —0.5 GeV^. Since the parton densities are defined at rather low momentum scale Qq = 
0.23 GeV^ we evolve them up to the experimental scale using the formalism of [Q. We concentrate 
mostly on the cross sections proportional to the imaginary part of the DVCS amplitude. Since Ei 
is concentrated in the ER-BL region with zero at \t\ = ^, Im£^i vanishes at all order. 

Similarly to the previous section we calculate the charge asymmetry for unpolarized exper- 
iments given in Eq. (ISOf ). As we see form the Fig. |^(a) it reaches the level of 5 — 15% for 
-A2 = 0.5 GeVl The single lepton spin asymmetry (|60|) is much larger and can be as big as 
20 — 30% [see Fig. ^(b)] which gives promises to measure it experimentally. Note that in both 



18 



cases the cross sections are very sensitive to the hehcity flip distribution E which is responsible 
for the orbital momentum of partons in the proton. 

In view of limited space let us consider finally only the spin asymmetry for a longitudinally 
and transversally polarized proton (0). In the former case azimuthal averaging of Eq. ( ^T]) is the 
same as in Eq. ( p(]| ) and leads to a sizeable asymmetry of order 20%. In the later case in order 
to extract the combination of SPDs multiplied by cos (0^/2) we define the following azimuthal 
asymmetry 



_ lit d<Pr^ - d<Pr^ ^^^^ 



The numerical estimate presented in Fig. ^ (d) demonstrates that it has a sizable effect which in 
contrast to the other symmetries does not vanish at the kinematical boundery = A ^j„ . 



5 Conclusions. 

In this paper we have given theoretical predictions for diverse asymmetries which can be measured 
in exclusive leptoproduction experiments of a real photon. Our estimates are rather encouraging, 
since they demonstrate a possibility to separate the contributions coming from different leading 
twist-two DVCS amplitudes by means of polarized lepton beam, and longitudinally and transver- 
sally polarized targets. This is the first step on the way to constrain the form of the SPDs 
from experimental data. The models used for numerical estimates lead to large charge and single 
(lepton) spin asymmetries of order 20%. However, since an asymmetry gives information on a 
linear combination of SPD one has to resort to other combinations of cross sections in order to 
disentangle a given distribution. 

We have not discussed in the present paper the NLO correction to the DVCS amplitudes. 
However, as has been shown in Ref. the latter could be very sizable and therefore change the 
LO predictions significantly in the valence quark region. Yet another important issue is the study 
of kinematical and dynamical higher twist corrections. Both of these problems deserve a detailed 
investigation and will be considered elsewhere. 

We would like to thank M. Amarian, M. Diehl, and A.V. Radyushkin for constuctive cor- 
respondence. This work was supported by DFG and BMBF. D.M. thanks G. Sterman for the 
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